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Description 

BACKGROUND OF THE INVENTION 

5 [0001] The present invention relates to the determination of properties of an optical device under test, e.g. the de- 
termination of the group delay of the optical device. The group delay is a fundamental property of optical devices, such 
as single mode optical fibers or optical components such as Bragg gratings which devices are used in the optical 
transmission of information. 

[0002] As for the background of optical transmission of data it has to be said that the premier feature of optical fiber 
10 js its extremely low loss. This has made it the dominant transmission medium for long link lengths. The loss charac- 
teristics of fiber determines were optical communication is practical. At 1 550 nanometer single mode optical fiber has 
an attenuation of 0.2 dB/km. This allows fiber optic signals to be propagated through very long length of fiber without 
regeneration. Telecommunication systems use the 1300 and 1550 nm windows for lowest loss in the fiber. Since a 
telecommunication system must cover a very large distance, the aforementioned attenuation of the single strength in 
'5 the fiber is of high importance. Therefore, the loss characteristics of optical fiber often limit the distance that a signal 
can propagate in the fiber. 

[0003] However, this is not always the case. In single mode fibers, chromatic dispersion can limit the distance over 
which fiber optic signals can propagate. Chromatic dispersion describes the fact that the speed of signal propagation 
in the fiber depends on the wavelength of the light. The consequence is that as the signal propagates through a long 

20 length of fiber, the edges of the waveform of the signal start to become more rounded. Eventually, the adjacent bits 
start to overlap in time causing the digital waveform to have poor readability. The amount of signal rounding depends 
beyond other parameters on the amount of chromatic dispersion in the cable. The problem of this pulse spreading is 
a problem in today's fiber communication because the increasing bit rate of the state of the art fiber communication 
systems, which bit rate reaches numbers of up to 40 Gbit/s per channel, brings the chromatic dispersion in the range 

25 of the bit resolution of such a 40 Gbit/s communication system. 

[0004] To compensate chromatic dispersion in a fiber chirped Bragg gratings were developed. In these gratings 
different wavelengths of the pulse need different amounts of time traveling through the grating to compensate the 
dispersion caused by the fiber. The concept of this compensation is schematically illustrated in Fig. 1 . Fig. 1 shows 
the original pulse 1 that is broadened to a spreaded pulse 2 while traveling through a fiber 3. With the help of a circulator 

30 4 a chirped Bragg grating 5 is introduced in the path of the pulse 2. As described before the resulting pulse 6 has the 
original shape of pulse 1 again. 

[0005] However, to use such a component as a chirped Bragg grating it is necessary to have exact knowledge of 
the chromatic dispersion caused by such a grating. The measurement of chromatic dispersion is accomplished by 
analyzing the group delay through the fiber as a function of wavelength. To measure the group delay a wavelength 

35 tunable optical source such as a tunable laser is used to generate coherent light at different wavelengths. In the prior 
art the wavelength of the tunable laser is then incremented step by step and for each wavelength step a group delay 
is determined. Finally, the group delay details are used to calculate the chromatic dispersion coefficient. The disad- 
vantage of measuring the group delays step by step is the time, which is needed to perform all the wavelength steps 
and the respective measurements. In other words this means it takes quite a long time to get data, which are so precise 

40 that this precision is sufficient for the needs of the telecommunication industry. 

SUMMARY OF THE INVENTION 

[0006] Therefore, it is an object of the invention to provide an improved determination of a property of an optical 
45 device under test, e.g. the group delay of the optical device. The object is solved by the independent claims. 

[0007] The main advantage of the present invention is that it is possible to measure the group delay and therefore 
the chromatic dispersion of an optical device, such as a fiber or a chirped Bragg grating, in a short amount of time 
while still keeping the precision of the data high enough for the needs of the telecommunication industry. This goal is 
reached by the invention by tuning the frequency of the coherent light beam of the laser from a maximum to a minimum 
50 of a given frequency range in a given time interval, e.g. making a wavelength sweep through the given wavelength 
range. Because of the inventive correction of effects caused by a non-linearity in the tuning of the laser this sweep can 
be done without detrimental broadening of the resulting spectral width. Method and apparatus of the invention avoid 
the aforementioned problems of the prior art and provide for exact data well keeping the measurement time low. 
[0008] In a preferred embodiment of the Invention there is performed a filtering of the Fourier transformed first signal 
55 with a high pass filter. It is further preferred to use a Hanning window as a shape for the high pass filter. Using this filter 
makes sure that a good elimination of not usable data is possible while still having important parts of data at the edges 
of the spectrum within the filter. Moreover, it is preferred to use a half of a Hanning window as a shape for the high- 
pass filter. 
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[0009] In a further preferred embodiment of the invention this high pass filter is adapted to the precision and the 
shape of the resulting spectrum of the signal corrected for a non-linearity of the laser by making an interferometric 
signal out of the corrected first phase signal, Fourier transforming the interferometric signal to get a spectral signal, 
determining a fraction of the maximum of the spectral signal, determining the abscissas of the intersections of the 

5 ordinate of the fraction which the curve of the spectral signal, determining the mean frequency f mean as the average 
of the abscissas, band pass filtering the spectral signal with a band pass filter having its center at the mean frequency 
and having a width greater than the width of the frequency range. The inventive concept of adaptive filtering of the 
signal corrected for a non-linearity of the laser is doubling the inventive success since by eliminating the non-linearity 
the resulting peak of the Fourier transformed spectrum is very sharp so that it can be used a filter width which is very 

io much smaller than known filter width in the prior art. As can be seen in the following detailed description of the invention 
and the respective drawings this inventive concept is resulting in very smooth curves showing the group delay. 
[0010] In a preferred embodiment the inventive method comprises the following steps: 

- splitting the first initial light beam into a first light beam and a second light beam, 

15 

- coupling the first light beam into the optical device under test, 

letting the second light beam travel a different path as the first light beam, 

20 - superimposing the first and the second light beam to produce interference between the first light beam and the 
second light beam in a resulting first superimposed light beam, 

detecting as a first signal the power of the first superimposed light beam as a function of time when tuning the 
frequency of the coherent light beam from a minimum to a maximum of a given frequency range in a given time 
25 interval, 

- splitting the second initial light beam in a third light beam and a fourth light beam, 

superimposing the third light beam and the fourth light beam after each light beam has traveled a different path, 
30 to produce interference between the third and the fourth light beam in a resulting second superimposed light beam, 

detecting as a second signal the power of the resulting second superimposed light beam as a function of time 
when tuning the frequency of the coherent light beam from a maximum to a minimum of a given frequency range 
in a given time interval, 

35 

using the detected second signal for deriving a non-linearity information about a non-linearity in a tuning gradient 
of the frequency when tuning the frequency of the coherent light beam from the maximum to the minimum of the 
given frequency range, and 

40 - using the non-linearity information for correcting effects on the first signal caused by the non-linearity to get a 
corrected first signal. 

[001 1 ] It is further preferred to flatten the first and/or the second signal before transf orming it by approximating it with 
polynoms of at least second order. 
45 [0012] Additionally it is preferred to tune the frequency with a tuning velocity of about 5 nm/s. 

[0013] Additionally it is preferred to tune the frequency within a tuning interval of about c/1400nm to c/1 600nm, c = 
velocity of light. 

[0014] Additionally it is preferred to detect the power at a rate of about 1 00kHz. 

[0015] The term "coherent" in this application means that the coherence length of the light beam is bigger than the 
so difference of the path of the first and second and third and fourth light beams, respectively. 
[0016] Other preferred embodiments are shown by the dependent claims. 

[0017] It is clear that the invention can be partly or entirely embodied or supported by one or more suitable software 
programs, which can be stored on or otherwise provided by any kind of data carrier, and which might be executed in 
or by any suitable data processing unit 

55 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0018] Other objects and many of the attendant advantages of the present invention will be readily appreciated and 
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become better understood by reference to the following detailed description when considering in connection with the 
accompanied drawings. The components in the drawings are not necessarily to scale, emphasis instead being placed 
upon clearly illustrating the principles of the present invention. Features that are substantially or functionally equal or 
similar will be referred to with the same reference sign(s). 

5 

Fig. 1 shows a chirped Bragg grating in a light path, 
Fig. 2 shows an embodiment of the invention, 
10 Fig. 3 shows a flow chart of the steps of an embodiment of the inventive method, 
Fig. 4 shows an illustration of the step 1 00 of Fig. 3, 
Fig. 5 shows an illustration of the step 1 06 of Fig.3, 

15 

Fig. 6 shows an illustration of the step 11 6 of Fig.3, 
Fig. 7 shows two possible shapes for filters used in step 11 6 of Fig. 3, 
20 Fig. 8 shows an illustration of the step 1 00 of Fig. 3, and 

Fig. 9 shows the group delay of the grating 1 2 of the example of Fig. 1 . 
DETAILED DESCRIPTION OF THE INVENTION 

25 

[001 9] Referring now in greater detail to the drawings, Fig. 2 shows a schematic illustration of a preferred embodiment 
of an apparatus 1 0 for interferometric determination of the group delay of a chirped Bragg grating 1 2 as a device under 
test. However, with the present invention it is possible to determine optical properties, e.g. the group delay ,of other 
devices, e.g. fibers or even air, also. 

30 [0020] The apparatus 1 0 comprises a tu nable laser 1 4 as a signal source for a coherent laser beam 1 6 , wherein the 
laser 14 is tunable between 1400 nm and at most 1600nm. The tunable laser 14 can be tuned continuously within a 
given frequency range in a given time interval. The laser beam 16 is coupled into a first beam splitter 18, which splits 
the coherent laser beam 16 in a first initial beam 20 and a second initial beam 22. The first initial beam 20 is coupled 
into a second beam splitter 24. With the help of the second beam splitter 24 the first initial beam 20 is split into a first 

35 beam 26 and a second beam 28. The first beam 26 is traveling a different and longer path than the second beam 28, 
since the first beam 26 is coupled via a circulator 52 into the grating 12. Reflected by the grating 12 the first beam 26 
is coupled via the circulator 52 in the first path 32 of the first beam 26 again. Then the first beam 26 and the second 
beam 28 are united with a third beam splitter 34 by superimposing the first beam 26 and the second beam 28 to a first 
resulting beam 36. The first resulting beam 36 is coupled into a first photodiode 38. The photodiode 38 transmits its 

40 output as a first signal to an analog/digital-converter (ADC) 40 (which can be National Instruments AD-MIO-1 6DE-1 0). 
The ADC 40 is connected to a transient recorder 42 that can be a FAST TR1202 of Fast Comtec GmbH, Germany. 
The laser 1 4 can be a H P81 680A of the applicant. The transient recorder 42 is coupled to a computer 44 as an evaluation 
unit of apparatus 10. 

[0021] Second beam splitter 24, first beam 26, second beam 28 and third beam splitter 34 built up a first Mach- 

45 Zehnder-interferometer. 

[0022] Exiting the first beam splitter 1 8 is a second initial beam 22. A fourth beam splitter 46 splits the second initial 
beam 22 into a third beam 48 traveling a third path and a fourth beam 50. The third beam 48 is traveling a different 
and longer path than the fourth beam 50, wherein the path length of the third path differs from the path length of the 
fourth path by more than 1 m, preferred by less than 20 m, more preferred by less than 9 m, even more preferred 

so between 7 and 9 m. This different path is symbolized by some loops 30. The third 48 and the fourth 50 beam are 
coupled into a fifth beam splitter 54 which units the third beam 48 and the fourth beam 50 to a superimposed second 
resulting beam 56. The second resulting beam 56 is then coupled into a second photodiode 58. The second photodiode 
58 transmits its output as a second signal to an analog/digital converter ADC 60, similar to the ADC 40. The ADC 60 
is also connected to the transient receiver 42. Fourth beam splitter 46, third beam 48, fourth beam 50 and fifth beam 

55 splitter 54 built up a second Mach-Zehnder interferometer. 

[0023] The system for determination of the group delay of the optical component 1 2 according to the present invention 
is based on the use of aforementioned Mach-Zehnder interferometers, both of which are supplied with light of the same 
light source 14. It is necessary to make the determination by interferometric measurements since the phase charac- 
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teristics of the component 12 can not be detected directly. 

[0024] The calculation of the properties based on the detected first and second signals is as follows: 

[0025] The laser 1 4 emits light 1 6 with field strength E dependent on the time t according to the following formula: 



E w (0=E 0 exp(y In \v mon {f)df) (3.1) 



10 with 



v_ 2ndt 



(3.2) 



*5 [0026] into the Mach-Zehnder interferometers, x> mom denoting the instantaneous frequency of the laser 14 . 

[0027] In the embodiment of Fig. 2 the frequency of the laser 1 4 is changed in a linear way so that the instantaneous 
frequency is 



20 



v marn(0=V 0 +Y-r 



(3.3) 



[0028] To calculate as a second signal the power of the resulting signal 56 of the second Mach-Zehnder interferometer 
detected with the second detector 58, one puts for the detected voltage U det that the light of the longer path is delayed 
by t = AL n/c when arriving at the fifth beam splitter 54. n is the refractive Index of the medium, AL is the difference of 
25 the path length and c is the velocity of light. 

[0029] The detected voltage is determined by the two superimposed fields that travel through the different paths of 
the interferometer: 



30 



40 



=2 
= 2 



exp(y -2* Jv^/VO+expO- In jv mom (t')df) 
o o 

l+cos(2;r jv moa {t<)dt>-2n \v mom (t')dt') 

o 0 
* + r "I 

1+008(2* jv n<m (t')dt')\ 
[0030] Only regarding of the changing part of this signal and putting r =fo(t) dt the result is: 



(3.4) 



45 



C/i, (t)=2 cos[2K (T^it + r)-r_(0)]. (3.5) 



[0031] Assuming t«t and approximating with differential quotient one reaches: 



50 
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10 



^(0=2-cos[2. r - (/ + r r ) - r - (0 r 
=2cos[2/r- ^Jt) t] (3.7) 

y ^average over r 



(3.6) 



[0032] The second signal of the second interferometer is therefore ideally, i.e. with a linear rise in frequency of the 
laser 14 over time a cosine, the frequency of which is determined by the length AL and the sweep velocity of the 
wavelength over time. 

[0033] Because the laser 14 will produce not an ideal sweep over the given frequency range and will produce non 
is linearities in this sweep one can expect a cosine which has a frequency which is changing over time. With this second 
signal it is possible according to the present invention to correct the first signal of the first interferometer for this non- 
linearities in therefore getting a resulting corrected first signal as if the laser 1 4 had performed a linear sweep. 
[0034] In the first interferometer in which the grating 1 2 is incorporated the group delay of this grating 1 2 is determined. 
For this purpose as the first signal the detected power of the first photodiode 38 is calculated in an analog way as 
20 explained above. However, in the following the first signal is only determined in its dependency of the optical frequency 
x>, where as the time dependency is determined later on. 

[0035] For the detected signal at the photodiode 38 the following formula applies: 

25 CdetM -lexp(./p DU7 <v).L Dar ) + exp(-/.p(v).L)l 2 =2[l4COs(P DUT <v)f. D ^. p(v)L)] (3.8) 

[0036] Whereas p is the phase constant of the second path, whereas p DUT is the phase constant of the first path 
through the grating 1 2. In this formula the fact is neglected that the grating 1 2 can have a reflection coefficient smaller 
than one. However, this has no importance for the changing part of the detected signal 36. For this changing part it 
30 can be written with <|>=p-L: 

^,(» / )=2 cos(^ iy7 (v)-^(v)) (3.9) 

35 

[0037] Hilbert-transforming this signal gives the analytic signal: 



(y)=i <*v(<i>DVT (y)-<fi(y)) <3.io> 

and the phase information: 

W V >=W V H( V ) < 3 - 11 > 

[0038] To determine the time dependency of the signal when performing a frequency sweep of the laser 14 and to 
find an evaluation rule for the group delay T g DUT , the functions § and fy^^y are now examined. 
[0039] For the phase of the path of the beam 28 without the grating 1 2 which path is assumed to be dispersion free 
it is possible to perform a Taylor approximation of the first order: 



55 
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(v) =^(v 0 )+^v>=v 0 (v-i/ 0 ) 
dv 

=f(v 0 )+2x^\a>=a> 0 -{v-v 0 ) (3/j2) 
=<*(v 0 ) + K -^-2n(v-v 0 ) 



where L is the length of the path without the grating 12. 

[0040] The same goes for the path 32 of the beam 26 with the grating 1 2: 



four (") =<W('o)+-y i2L |»' = V(»'-Vo) 

1^ (313) 
+ l^rL = l , o . (v _ Vo) 2 +% 

2 */v 1 

where y is denoting the residual error. Performing the same steps as in formula 3.1 2 one reaches the following formula: 



+ \i^ {Tdvt (v)) I v=v ° - 2ni y- v ^ +4/ (3 ' 14> 



with 



dr pur = dx Dm dX = dv DUT 1 _ dz X 2 
dv dX dv dX dX c' " 

dX{x) 



[0041] By putting this formula into 3.11 this gives for the phase of the analytic signal: 



+ 2n[z D(/r (v 0 )-z(v 0 )]-(v-v 0 ) 

dr.. X 2 . (3-16) 

-— i^o — *(v-Vo) 
dX c 



= A<|> 0 (3.17) 



+2kAt 0 (v-v 0 ) 



(3.18) 
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(3.20) 

[0042] To determine the group delay of the grating 12 it is sufficient to differentiate the face <|> ctet with respect to <o. 
This gives: 



10 



15 



^T-] (3.21) 



dco In dv 2/r dv dv 



[0043] This is the information searched for since only the behavior of the group delay with respect to the frequency 
is interesting, rather than an absolute value. 

[0044] However, it is easy to extract from formulas 3.1 7 - 3.20 information about a deviation of this linear rise of the 
20 group delay with respect to the wavelength for gratings with linear chirp. This goes by approximating the measured 
phase signal with a polynomial of second order and subtracting this polynomial of the phase signal. After this only the 
part 3.20 remains which represents the deviation of the group delay with respect a linear chirp grating. 
[0045] With the help of the diagram of Fig. 3 the main steps of an embodiment of the inventive method are described, 
whereas some more details of these steps follow later on: 
25 [0046] In Fig. 3 ovals describe signals and boxes describe measures of the evaluation unit 44. First of all the second 
or reference signal 56 and the first signal 36 are towed to the zero line in step 1 00 (see Fig. 4). The resulting signals 
102 and 104 are now oscillating about the zero line. In step 106 any non-linearities in the sweep of the laser 14 are 
corrected with the help of the second signal 56 (see Fig. 5). Resulting is a signal 108 corrected for any non-linearities 
in the laser sweep. With the help of the aforementioned formulas are evaluated in step 1 1 0 the group delay 112 and 
30 the linear deviation of linear chirp 114. 

[0047] With the help of a filter step 116 it is possible to get filtered signals 1 1 8 and 1 20 of the group delay and of the 
deviation from the linear chirp of the grating 12. 

[0048] In Fig. 4 the principle of the function 1 00 is graphically illustrated. A is the amplitude of the signals 36 and 56. 
152 denotes the reflection spectrum of the grating 12. 124 denotes the ideal center line of the spectrum 122. 
35 [0049] The function 1 00 is illustrated in Fig. 8 and comprises the following steps: 

making the signal 36, 56 oscillating about a zero line by: 

determining the points of mean value of the signal 36, 56 by: 

40 

determining the maximum and the minimum of the signal 36, 56 in a predetermined first range of time 
smaller than the total range of time, 

- determining a mean value between the maximum and the minimum, 

determining the maximum and the minimum of the signal 36, 56 in a predetermined next range of time 
adjacent the already examined range of time, 

- determining a mean value between the maximum and the minimum, 
repeating the last two steps until the complete time interval is covered, 

choosing the predetermined range of time by: 
55 - determining the average period of the oscillations of the signal 36, 56, 

choosing the size of the range so that more than two average periods fit in the chosen range of time. 



45 
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[0050] Figure 5 illustrates a part of the function of step 106 to correct unwitting non-linearities in the sweep of the 
laser 14. The step 106 comprises the following steps (in the following the symbol T is used instead of V"): 

- deriving the non-linearity information by: 

5 

- transforming the second signal 56 to get a Fourier transformed second signal, 

eliminating the negative parts of the Fourier transformed second signal to get a non-negative Fourier trans- 
formed second signal, 

10 

retransforming the non-negative Fourier transformed second signal to get an analytic signal of the second 
signal 56, 

- determining the phase of the analytic signal to get as a second phase signal the phase as a function of time 
15 of the second signal 56, 

using the second phase signal for determining as the non-linearity information the frequency as a function of 
time of the second signal 56, 

20 - deriving a first phase signal by: 

- transforming the first signal 36 to get a Fourier transformed first signal, 

eliminating the negative parts of the Fourier transformed first signal to get a non-negative Fourier transformed 
25 first signal, 

retransforming the non-negative Fourier transformed first signal to get an analytic signal of the first signal 36, 
determining the phase of the analytic signal to get as a first phase signal as a function of time of the first signal, 

30 

- correcting the effects on the first signal 36 caused by the non-linearity by: 

using the non-linearity information to interpolate the first phase signal of the first signal 36 on a linear scale of 
frequencies to get a corrected first phase signal, 

35 

[0051] For determining the non-linearity information the following steps are performed: 

- determining the frequency f(n) of the second signal 56 as a function of n discrete points of time, n= 1 , N, on the 
basis of the second phase signal by: 

40 

determining the second phase signal <|>(n) at the n points of time, 

determining the maximum 4> max of the second phase signal, 

45 - using a predetermined maximum frequency f max of the frequency range, a predetermined average tuning 

velocity during tuning the frequency and the maximum <|> max of the second phase signal to determine for each 
of the n points of time the frequency f(n) according to the formula: f(n)-[(f max -f mln )/(fr max ]<Mn) ( 

[0052] For the interpolating step, the following sub-steps are performed: 

50 

- getting the linear scale f Ijn (n) of frequencies by: 

- using the predetermined maximum frequency f^ of the frequency range and the predetermined minimum 
frequency f min of the frequency range to determine the linear scale f lin (n) of frequencies according to the for- 

55 mula:f lhl (n)=[(f max -f mln )/(N-1)]n, 

preferably sorting the absolute values of f (n) monotonically, 



9 



EP1 191 321 A1 



using f (n) for interpolating the first phase signal of the first signal 36 on the linear scale of frequencies f Nn (n). 
The interpolation can be done linearly or with the help of polynoms or splines. 

[0053] The resulting signal 1 08 which Is corrected for any non linearities in the sweep of the laser 1 4 is used in step 
5 110 for evaluate the group delay 112 and the deviation 114 from the linear chirp of the grating 12. This evaluation is 
done with the help of the above-mentioned mathematical formulas and can comprises deriving a group delay 112 of 
the optical device under test 12 by differentiating the compensated first phase signal with respect to the frequency. 
[0054] The method can further comprise the following steps: 

10 - Ignoring at the begin of the tuning a predetermined amount of values of the corrected first phase signal to eliminate 
teething troubles out of the corrected first signal 1 08, 

approximating the group delay 112 with polynoms of at least second order to get an approximated group delay, 

15 - subtracting the approximated group delay from the group delay 1 1 2 to get a non-linear part of the group delay 1 1 2. 

using the non-linear parts of the group delay 1 1 2 to determine the mean signal power of a deviation from a linear 
of the group delay 112 of the device under test 12. 

20 - using the square coefficient of the polynomial to determine the mean gradient of the group delay 112. 

[0055] Having evaluated the group delay 112 and the deviation 114 it is possible but not necessary to make a filter 
step 1 1 6 to filter the signals 112,114. Actually, this filter step is preferably performed before the calculating of the group 
delay 112, as follows: 

25 

adapting the filtering to the shape of the corrected first signal 1 08 by: 

a: making an interferometric signal out of the corrected first phase signal, 
30 - b: Fourier transforming the interferometric signal to get a spectral signal, 

- c: determining a fraction, preferred the half, of the maximum of the spectral signal, 

d: determining the abscissas of the intersections of the ordinate of the fraction with the curve of the spectral 
35 signal, 

- e: determining the mean frequency f mean as the average of the abscissas, 

- f: band-pass filtering the spectral signal with a band-pass filter having its center at the mean frequency and 
40 having a width greater than the width of the frequency range, 

determining the width of the band-pass filter by: 

- a: predetermining or estimating the maximal range GD^^ of the group delay, 

45 

- b: determining the mean value GD mean of the group delay according to the steps c-e of claim 26, 

- c: calculating the filter width according to the formula: filter width= 1 mean (GO mngB /GD meai J, 
so - filtering with a filter having the shape of a Hanning window. 

[0056] Resulting is a filtered group delay 118 and a filtered deviation 120. 

[0057] Figure 6 shows the principle of determination of the center frequency f mean of the interferometric spectrum. 
[0058] Figure 7 shows a graph showing different filter shapes for use with the inventive method, a flat top filter shape 
55 130 and a Hanning form filter shape 140. Figure 9 shows the group delay 112 over wave length of a chirped Bragg 
grating 1 2, which is a FBGDC-HP1 A of the company JDS Fitel. It is clearly seen that the graph of the group delay 1 1 2 
is very smooth with respect to graphs determined with methods of the prior art. 
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Claims 

1 . A method of determination of a property of an optical device under test (12), comprising the steps of: 
using a first initial coherent light beam (20), 
changing a first initial property of the first initial light beam (20), 
coupling the first initial light beam (20) to the device under test (12), 

detecting a first signal of the first initial light beam (20) received from the device under test (12), 

correcting a non-linearity in the first signal caused by a non-linearity in the change of the first initial property 
by interpolating the first signal on a linear scale. 



2. The method of claim 1 , further comprising the steps of: 

using a second initial coherent light beam (22), 

20 - changing a second initial property of the second initial light beam (22), 

detecting a second signal of the second initial light beam (22) without coupling it to the device under test (1 2), 
to discover a non-linearity in the second signal caused by a non-linearity in the change of the second initial 
property, 



using the discovered non-linearity of the detected second signal to interpolate the first signal. 

3. The method of any one of the claims 1-2, further comprising the steps of: 
30 - producing a coherent light beam (1 6), 

- splitting the coherent light beam (1 6) into a first initial light beam (20) and a second initial light beam (22). 

4. The method of any one of the claims 1 -3, further comprising the step of: 

35 

detecting the first resulting property simultaneously with the second resulting property. 

5. The method of any one of the claims 1 -4, 

40 - changing the first initial property simultaneously with the second initial property. 

6. The method of any one of the claims 1 -5, 

wherein the first initial property and the second initial property being the same initial property. 

45 7. The method of any one of the claims 1-6, 

wherein the initial property being the frequency of the coherent light beam (16). 

8. The method of any one of the claims 1 -7, further comprising the steps of: 

so - transforming the first signal in a number of phase signals over a linear scale of a number of points of time, 

- transforming the second signal in a number of frequency signals over the same linear scale of points of time 
to discover a non-linearity in the second signal caused by a non-linearity in the change of the initial property, 
the initial property being the frequency of the coherent light, 



55 



assigning the transformed first signal to the transformed second signal, 
interpolating the assigned transformed first signal on a linear scale of frequencies. 
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9. The method of claim 8 f further comprising the steps of: 

- creating the linear scale of frequencies f„ n (n) according to the formula f Dn (n) = (f min - f max ) x (n/N), n Element 
1 N, N being the number of points of time. 

5 

10. The method of any one of the claims 1 -9, further comprising the steps of: 

- splitting the first initial light beam (20) into a first light beam (26) and a second light beam (28), 
10 - coupling the first light beam (26) to the optical device under test (1 2), 

- letting the second light beam (28) travel a different path as the first light beam (26), 

superimposing the first (26) and the second (28) light beam to produce interference between the first light 
15 beam (26) and the second light beam (28) in a resulting first superimposed light beam (36), 

detecting as a first signal the power of the first superimposed light beam (36) as a function of time when tuning 
the frequency of the coherent light beam (1 6) from a minimum to a maximum of a given frequency range in a 
given time interval, 



20 



55 



splitting the second initial light (22) beam in a third light beam (48) and a fourth light beam (50), 



- superimposing the third light beam (48) and the fourth light beam (50) after each light beam (48, 50) has 
traveled a different path, to produce interference between the third (48) and the fourth (50) light beam in a 
25 resulting second superimposed light beam (56), 

detecting as a second signal the power of the resulting second superimposed light beam (56) as a function of 
time when tuning the frequency of the coherent light beam (16) from a maximum to a minimum of a given 
frequency range in a given time interval, 

30 

using the detected second signal for deriving a non-linearity information about a non-linearity in a tuning gra- 
dient of the frequency when tuning the frequency of the coherent light beam (16) from the maximum to the 
minimum of the given frequency range, and 

35 - using the non-linearity information for correcting effects on the first signal caused by the non-linearity to get a 

corrected first signal (1 08). 

1 1 . The method of any one of the claims 1 -1 0, further comprising the steps of: 

40 - deriving the non-linearity information by: 

- transforming the second signal to get a Fourier transformed second signal, 

eliminating the negative parts of the Fourier transformed second signal to get a non-negative Fourier 
45 transformed second signal, 

retransforming the non-negative Fourier transformed second signal to get an analytic signal of the second 
signal, 

so - determining the phase of the analytic signal to get as a second phase signal the phase as a function of 

time of the second signal, 

using the second phase signal for determining as the non-linearity information the frequency as a function 
of time of the second signal. 



12. The method of any one of the claims 1-11 , further comprising the steps of: 
deriving a first phase signal by: 
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- transforming the first signal to get a Fourier transformed first signal, 

eliminating the negative parts of the Fourier transformed first signal to get a non-negative Fourier trans- 
formed first signal, 

retransforming the non-negative Fourier transformed first signal to get an analytic signal of the first signal, 

determining the phase of the analytic signal to get as a first phase signal as a function of time of the first 
signal. 

13. The method of any one of the claims 1-12, further comprising the steps of: 

correcting the effects on the first signal caused by the non-linearity by: 

using the non-linearity information to interpolate the first phase signal of the first signal on a linear scale 
of frequencies to get a corrected first phase signal. 

14. The method of any one of the claims 1 -13, further comprising the steps of: 

determining the frequency f(n) of the second signal as a function of n discrete points of time, n= 1, N, on 
the basis of the second phase signal to determine the non-linearity information by: 

determining the second phase signal <|>(n) at the n points of time, 

determining the maximum <fr max of the second phase signal, 

using a predetermined maximum frequency f max of the frequency range, a predetermined average tuning 
velocity during tuning the frequency and the maximum of the second phase signal to determine for 
each of the n points of time the frequency f(n) according to the formula: ^nHAnax^iriJ^muMt 11 )- 

15. The method of any one of the claims 1 -14, further comprising the steps of: 

getting the linear scale fnn(n) of frequencies by: 

using the predetermined maximum frequency f max of the frequency range and the predetermined minimum 
frequency f mln of the frequency range to determine the linear scale f Hn (n) of frequencies according to the 
formula: WnHlUWfN-l )]n, 

preferably sorting the absolute values of f (n) monotonically. 

16. The method of any one of the claims 1 -15, further comprising the steps of: 

using f(n) for interpolating the first phase signal of the first signal on the linear scale of frequencies f, in (n). 

17. The method of any one of the claims 1-16, further comprising the steps of: 

deriving transmissive and/or reflective properties of the optical device under test from the compensated first 
signal. 

18. The method of any one of the claims 1 -17, further comprising at least one of the following steps of: 

- deriving a group delay (112) of the optical device under test (12) as a function of frequency from the corrected 
first signal (108), 

deriving the chromatic dispersion coefficient of the optical device under test (12) as a function of frequency 
from the corrected first signal (108). 

19. The method of any one of the claims 1 -1 8, further comprising the steps of: 
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- deriving a group delay (112) of the optical device under test (12) by differentiating the corrected first phase 
signal with respect to the frequency. 

20. The method of any one of the claims 1-19, further comprising the steps of: 

5 

ignoring at the begin of the tuning a predetermined amount of values of the corrected first phase signal to 
eliminate teething troubles out of the corrected first signal (1 08). 

21 . The method of any one of the claims 1 -20, further comprising the steps of: 

10 

- approximating the group delay (112) with polynoms of at least second order to get an approximated group delay, 

subtracting the approximated group delay from the group delay (112) to get a non-linear part of the group 
delay (112). 

15 

22. The method of any one of the claims 1 -21 , further comprising the steps of: 

using the non-linear part of the group delay to determine the mean signal power of a deviation from a linear 
group delay of the device under test (12). 

20 

23. The method of any one of the claims 1 -22, further comprising the steps of: 

using the square coefficient of the polynomial to determine the mean gradient of the group delay (112). 
25 24. The method of any one of the claims 1 -23, further comprising the steps of: 

- making the first signal oscillating about a zero line by: 

- determining the points of mean value of the first signal, 

30 

interpolating a curve through these points, 

- subtracting the values of the curve from the first signal to get a corrected first signal oscillating about the 
zero line. 

35 

25. The method of any one of the claims 1 -24, further comprising the steps of: 

- determining the points of mean value by extracting all points with a maximum gradient. 
40 26. The method of any one of the claims 1 -25, further comprising the steps of: 

making the first signal oscillating about a zero line by: 

- determining the points of mean value of the first signal by: 

45 

determining the maximum and the minimum of the first signal in a predetermined first range of time 
smaller than the total range of time, 

- determining a mean value between the maximum and the minimum, 

50 

determining the maximum and the minimum of the first signal in a predetermined next range of time 
adjacent the already examined range of time, 

determining a mean value between the maximum and the minimum, 

55 

repeating the last two steps until the complete time interval is covered. 
27. The method of any one of the claims 1 -26, further comprising the steps of: 
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choosing the predetermined range of time by: 

determining the average period of the oscillations of the first signal, 
5 - choosing the size of the range so that more than two average periods fit in the chosen range of time. 

28. The method of any one of the claims 1-27, further comprising the steps of : 
■ determining the points of mean value by: 

10 

determining the maximum of the Fourier transformed signal of the first signal, 
using the maximum to determine a size of a high-pass filter, 
is - filtering the Fourier transformed first signal with the high-pass filter. 

29. A method of determination of a property of an optical device under test (12), comprising the steps of: 

detecting a change of a signal with time, being the basis for deriving the property, 

20 

filtering the detected signal by: 

transforming the detected signal to get a Fourier transformed signal, 
25 - filtering the Fourier transformed signal with a filter to get a filtered Fourier transformed signal, 

retransforming the filtered Fourier transformed signal to get a filtered signal, 
deriving the property on the basis of the filtered signal. 

30 

30. The method of claim 29, further comprising the steps of: 

the detected signal being corrected for a non-linearity by the method of any one of the claims 1 - 28 to get a 
corrected first signal (108). 

35 

31. The method of any one of the claims 1 -30, further comprising the steps of: 

using the corrected first signal to calculate the corrected first phase signal versus frequency, 

40 - filtering the corrected first phase signal by: 

Hilbert transforming it before filtering it to get a corrected signal to be filtered according to the steps of 
any one of the claims 29 or 30. 

45 32. The method of any one of the claims 1 -31 , further comprising the steps of: 

- filtering the corrected first signal (1 08) before calculating the group delay (112). 
33. The method of any one of the claims 1 -32, further comprising the steps of: 

50 

- adapting the filtering to the shape of the corrected first signal (1 08) by: 

a: making an interferometric signal out of the corrected first phase signal, 
55 . b: Fourier transforming the interferometric signal to get a spectral signal, 

- c: determining a fraction, preferred the half, of the maximum of the spectral signal, 
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d: determining the abscissas of the intersections of the ordinate of the fraction with the curve of the spectral 



e: determining the mean frequency f mean as the average of the abscissas, 

5 

- f : band-pass filtering the spectral signal with a band-pass filter having its center at the mean frequency 
and having a width greater than the width of the frequency range. 

34. The method of any one of the claims 1-33, further comprising the steps of: 

10 

- determining the width of the band-pass filter by: 

a: predetermining or estimating the maximal range GD^ge of the group delay, 
15 b: determining the mean value GDmean of the group delay according to the steps c-e of claim 33, 

- c: calculating the filter width according to the formula: filter width= f me an( GD rang^ GD mean)- 

35. The method of any one of the claims 1 -34, further comprising the steps of: 

20 

- subtracting a gradient in the group delay from the group delay, 
predetermining the maximum range GD^^ of the group delay, 

25 - performing the steps b and c of claim 34, 

- calculating the group delay, 

- adding the subtracted gradient to the calculated group delay. 

30 

36. A software program or product, preferably stored on a data carrier, for executing the method of one of the claims 
1 to 35 when run on a data processing system such as a computer. 

37. An apparatus of determination of properties of an optical device under test, comprising: 

35 

a first beam splitter (18) in a path of a coherent light beam (1 6) for splitting the coherent light beam (1 6) into 
a first initial light beam (20) traveling a first initial path and into a second initial light beam (22) traveling a 
second initial path, 

40 - a second beam splitter (24) in that first initial path for splitting the first initial light beam (20) into a first light 

beam (26) and traveling a first path (32) and into a second light beam (28) traveling a second path, 

- a place in that first path (32) for coupling the first light beam (26) to the optical device under test (12), 

45 a third beam splitter (34) in that first (32) and in that second path for superimposing the first (26) and the 

second (28) light beam after the second light beam (28) has traveled a different path as the first light beam 
(26) to produce interference between the first light beam (26) and the second light beam (28) in a resulting 
first superimposed light beam (36) traveling a first resulting path, 

so a first power detector (38) for continuously detecting as a first signal the power of the first superimposed light 

beam (36) as a function of time when tuning the frequency of the coherent light beam (16) from a minimum 
to a maximum of a given frequency range in a given time interval, 

- a fourth beam splitter (46) for splitting the second initial beam (22) in a third light beam (48) traveling a third 
55 a path and a fourth light beam (50) traveling a fourth path, 

- a fifth beam splitter (54) in that third and in that fourth path for superimposing the third light beam (48) and the 
fourth light beam (50) after each light beam (48, 50) has traveled a different path, to produce interference 
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between the third (48) and the fourth light beam (50) in a resulting second superimposed light beam (56) 
traveling a second a resulting path, 

a second power detector (60) for continuously detecting as a second signal the power of the resulting second 
superimposed light beam (56) as a function of time when tuning the frequency of the coherent light beam (1 6) 
from a maximum to a minimum of a given frequency range in a given time interval, 

an evaluation unit (44) for deriving optical properties of the optical device under test (12), for using the detected 
second signal for deriving a non-linearity information about a non-linearity in a tuning gradient of the frequency 
when tuning the frequency of the coherent light beam (16) from the maximum to minimum of the given fre- 
quency range, and for using the non-linearity information for correcting effects on the first signal caused by 
the non-linearity to get a corrected first signal. 

38. The apparatus of claim 37, 

further having computing means capable of performing a least one of the further steps according to any of the 
claims 2-36. 

39. The apparatus of anyone of the claims 37 or 38, further comprising : 

a circulator (52) at that place in that first path (32) to enable the apparatus (1 0) to examine reflective optical 
components also. 

40. The apparatus of anyone of the claims 37-39, 

where the second beam splitter (24), the third beam splitter (34) and the first detector (38) build up a first Mach- 
Zehnder interferometer, 

and the fourth beam splitter (46) and the fifth beam splitter (54) and the second detector (58) build up a second 
Mach-Zehnder interferometer. 
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